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ABSTRACT 

We investigate the spectral and timing properties of the millihertz quasi- 
periodic oscillations (mHz QPOs) in neutron-star low-mass X-ray binary 4U 
1636-53 using XMM-Newton and Rossi X-ray Timing Explorer (RXTE) ob¬ 
servations. The mHz QPOs in the XMM-Newton/RXTE observations show 
significant frequency variation and disappear right before type I X-ray bursts. 
We find no significant correlation between the mHz QPO frequency and the 
temperature of the neutron-star surface, which is different from theoretical 
predictions. For the first time we observed the full lifetime of a mHz QPO 
lasting 19 ks. Besides, we also measure a frequency drift timescale ^15 ks, we 
speculate that this is the cooling timescale of a layer deeper than the burn¬ 
ing depth, possibly heated by the previous burst. Moreover, the analysis of 
all X-ray bursts in this source shows that all type I X-ray bursts associated 
with the mHz QPOs are short, bright and energetic, suggesting a potential 
connection between mHz QPOs and He-rich X-ray bursts. 


Key words: X-rays: binaries; stars: neutron; accretion, accretion discs; X- 
rays: individual: 4U 1636—53 


E-mail: m.lyu@astro.rug.nl 


2 Ming Lyu et al. 

1 INTRODUCTION 


The millihertz quasi-periodic oscillations (mHz QPOs) in neutron-star low- mass X-ray bina¬ 


ries ( 


MXBs) were first detected in 4U 1608-52, 4U 1636-53, and Aql X-1 by 


Revnivtsev et al 


fl200ll) . The mHz QPOs occur only within a narrow range of X-ray luminosities, -h 2 - 20 kev — 


(5 — 11) X 10^® ergs s in the frequency range 7 — 9 mHz with strong flux variation 


below 5 keV. Besides, t he mHz QPOs o 


flRevnivtsev et ah 


2001 


Altamirano et ah 


;en c 


200811 ■ 


isappear right before a type I X-ray burst 


It has been proposed that the mechanism responsible for the mHz QPOs is connected 
with s ome sort of nuclear burning on the neutron star surface. For instance. 


Revnivtsev et ah 


(120011) suggested that the mHz QPOs resulted from a special mode of nuclear burning on 
the neutron star su rface when the source was within a certain range of mass accretion rate. 

(1200711 further proposed that the mHz QPOs were due to marginally stable nu- 


Heger et al. 


clear burning of Helium on the surface of accreting neutron stars, with a predicted oscillation 
timescale of ~100 s, consistent with the observed ~2 minutes period of the mHz QPOs. 


Several results 


Yu fc van der Klis 


rave strengthened the nuclear burning explanation of the mHz QPOs. 


20021) detected mHz QPOs and kilohertz QPOs in the same observations 


of 4U 1608-52. The 2 — 5 keV X-ray count rate in the mHz QPO cycles was anti-correlated 
with the frequency of the kilohertz QPOs. This result was consistent with the nuclear burning 
origin of the mHz QPOs: As the luminosity increases, the stresses of radiation coming from 
the neutron star surface push the inner di sc outward so that the frequency of the kilohertz 


QPOs decreases. 


Altamirano et al. 


(120081 ) found a systematic frequency drift of the mHz 


QPO until the QPO disappeared before a type I X-ray burst in 4U 1636—53, indicating a 
connection between the mHz QPOs and the unstable nuclear reaction on the neutron star 
surface. 


Linares et al 


(l2010h reported a new detection of mHz QPOs at ~4.5 mHz in the neutron 
star transient IGR J17480-2446 at a persistent luminosity L 2 - 50 fceV ~ 10^® erg s“^, about 
one order of magnit ude higher than observed in the previous mHz QPO sources. In addition. 


Linares et al. 


(120121 ) found a transition between thermonuclear bursts and mHz QPOs in this 


source: the X-ray bursts gradually evolved into the mHz QPO as accretion rate increased, 
and vice versa. This behaviour was similar to the pr ediction frorn the o ne-zone models and 
simulation results of the marginally stable burning of lHeger et al.l (j2007ll . further conhrming 
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that the mHz QPO was connected with marginally stable bnrning on the nentron star 
snrface. 


In the model of lHeeer et al.l (120071) . the reqnired mass accretion rate for marginally stable 
bnrnin g is ~ Mf,hh, about ten times higher than the rate deduced by the ayerage X-ray lumi 


nosity (IRevnivtsev et ah 


2001 


Yu fc van der Klis 


2002 ; 


Altamirano et ah 


200811 ■ 


Keek et ah 


(120091) found that the critical accretion rate at which mHz QPOs were observed could be 
explained if they considered the turbulent chemical mixing of the accreted fuel together with 
a higher heat flux from the crust. Furthermore, the cooling of the deep l ayers also explains 


the observed frequency drifts before X-ray bursts. Recently, 


Keek et al, 


(1201411 investigated 


the influence of the chemical composition of the accreted fuel and nuclear reaction rates 
on the marginally stable burning process; their results showed that no allowed variation in 
accretion fuel composition and nuclear reaction rate was able to produce marginally stable 
burning at the observed accretion rate. 


Lyn et ahl (j2014al ) found that the frequency of the mHz QPOs correlates with the tem¬ 


perature of the neutron star surface in 4U 1636-53 using an XMM -Newton ob s ervati on. This 


frequency-temperature correlation agrees well with the results of 


Keek et al 


( 20091) : As the 


burning layers cool down, the 
until a burst happens. Besides, 


requency o 


the o scillations decreases by tens of per cents 


Lyn et al.l (l2014al ) found that the QPO in that observation 


disappeared after one burst, and it subsequently reappeared ~25.3 ks later, which was the 
longest reappearance time scale so far measured for a mHz QPO in any LMXBs. 


In this work we investigate the frequency-temperature correlation in the mHz QPO cycle 
using XMM-Newton observations of the LMXB 4U 1636-53 to take a closer look into the 
mechanism behind the mHz QPOs. Using simultaneous observations with the Rossi X-ray 
Timing Explorer (RXTE), this is the hrst time that we are able to carry out a wide energy 
band (0.8—100 keV) spectral analysis in the mHz QPO cycles. Furthermore, we explored 
the possible relation between the mHz QPOs and the associated X-ray bursts, which may 
provide a new way to investigate the origin of mHz QPOs. In Section 2 we describe the 
observations and data analysis, and in Section 3 we present our main results. Finally, we 
discuss the results in the framework of the marginally stable burning model in Section 4. 
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2 OBSERVATIONS AND DATA REDUCTION 


The three XMM-Newton observations of 4U 1636-53 that we use in this work w ere performed 
with the European Photon Imaging Camera, EPIC-PN fjStriider et ahll200lh in the timing 
mode, with a time resolution of 0.03 ms and in the energy band 0.5—12 keV. The three 
observations were taken between 2007 to 2009, with a total duration of about 90 ks. We 
selected RXTE observations that were taken simultaneo usly with those three XMM-Newton 


observations. The Proportional Counter A rray (PCA; 
Energy X-ray Timing Experiment (HEXTE; 


Jahoda et al. 


Rothschild et al. 


2006) and the High 


1998 1 onboard RXTE covered 


the energy range 2—60 keV and 15—250 keV, respectively. Both the PCA and HEXTE data 


were included in this work in order to get a wide energy 
the XMM-Newton observation (ObsID 0606070301) in 


pand cover a ge. We also re-analyzed 


Lyn et ahl (j2014al) with the latest 


calibration hies, and selected the simultaneous RXTE observation (ObsID 94310-01-03-000) 
for the spectral analysis. 


In this work, we used the Science Analysis System (SAS) version 14.0.0 for the XMM- 
Newton data reduction, with the latest calibration hies applied (as of April 2015). We selected 
the Rate-Dependent PHA (RDPHA) method in the command epproc for the correction of 
energy scale rate dependence in EPIC-PN Timing Mode exposures. We applied the task 
barycen to convert the arrival time of photons from the local satellite frame to the barycenter 
of the Solar system. We found that there was some moderate pile up in all XMM-Newton 
observations according to the results of the epatplot test. We hnally selected a 10-column 
wide region centered at the position of the source, excluding the central column for Obs 1 
and the three central columns for the other observations. We only selected single and double 
events (pattern ^ 4) to extract spectra and light curves. 


Furthermore, we investigated all X-ray bursts of 4U 1636-53 detected by the PCA 


flJahoda et al.M2006[l in order to explore the possible connection between mHz QPOs and 


the associated X-ray bursts. For this we produced 0.25-s light curves from the Standard-1 
da ta and searched for X-ray bursts in these light curves following the procedure described 


m 


Zhang et ahl (1201111 . A total of 338 type-I X-ray bursts have been detected in these data. 


We used the Standard-2 data to calculate the source X-ray colours. We dehned soft 
and hard X-ray colours as t he 3.5 — 6.0/2.0 — 3.5 keV and 9.7 — 16.0/6.0 — 9.7 keV count 


rate ratios, respectively (see 


Zhang et al 


2011 


, for details). We parametrized the position 


of the source on the colour-colour diagram (CCD) by the length of a solid curve (see. 
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Table 1. XMM-Newton / simutaneous RXTE observations of 4U 1636—53 used in this work. 


Observation 

Instrument 

Observation ID 

Start Time 

Observation Length (ks) 

Obs 1 

XMM-Newton 

0500350301 

2007-09-28 15:43 

30.6 


RXTE 

93091-01-01-000 

2007-09-28 14:48 

28.8 


RXTE 

93091-01-01-00 

2007-09-28 22:48 

9.5 

Obs 2 

XMM-Newton 

0500350401 

2008-02-27 04:14 

38.6 


RXTE 

93091-01-02-000 

2008-02-27 03:47 

28.8 


RXTE 

93091-01-02-00 

2008-02-27 11:47 

10.0 

Obs 3 

XMM- N e wton 

0606070101 

2009-03-14 15:11 

23.8 


RXTE 

94310-01-01-02 

2009-03-14 15:05 

2.4 


RXTE 

94310-01-01-04 

2009-03-14 16:40 

4.0 


RXTE 

94310-01-01-00 

2009-03-14 18:18 

13.4 

Obs 4* 

XMM- N ewton 

0606070301 

2009-09-05 01:57 

40.9 


RXTE 

94310-01-03-000 

2009-09-05 01:17 

23.2 


* We re-analysed the spectral properties of this observati on in this work. F or the timing properties of this observation, please 

refer to lLvu et al.l ll2014alV 



Figure 1. Colour-colour diagram of 4U 1636-53. Ea ch magenta point (grey in the black and white version of this figure) 
represents the averaged Crab-normalised colours fsee IZhang et aLll2Qll]. for details) of a single RXTE observation. In the 
legend we indicate the observations used in this work. The position of the source in the diagram is parameterised by the length 
of the blue solid curve S'a, the details of S'a curve is shown in the §2. 


e.g. 


Mendez et al. 


1999 


Zhang et al. 


201 ll) . fixing the values of S'a = 1 and Sa = 2 at the 


top-right and the bottom-left vertex of the CCD, respectively. 

The colour-colour diagram in Figure [U shows that the four XMM-Newton observations 
were performed when the source was in the transitional state between the hard and the soft 
spectral state. The details of the XMM-Newton and simultaneous RXTE observations in 
this work are listed in Table [1] 
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2.1 Timing data 


We extracted a 1-s resolution light curvqll in the 0.2—5 keV band and produced an av¬ 
erage power spectrum (Figure [2]) for each XMM-Newton observation, after excluding in¬ 


strum ent dropouts and X-ray bursts. We used the ISIS Version 1.6.2-27 flHouck fc Denicola 


2000l) to compute the Fourier transform to intervals of 512-s duration using the command 
sitar_avg_psd, in a frequency range of 1.95x10“^ to 0.5 Hz. To explore the frequency evo¬ 
lution of the mHz QPO, we generated a dynamic power spectrum with the frequenc y scale 


being oversampled by a factor of 100 using the Lomb-Scargle periodogram method (ILomb 


1976 


Scargle 


19821) . The dynamic power spectra of the four XMM-Newton observations are 


shown in Figure [3l each column in those plots represents the power spectrum extracted from 
a time interval of 1130 seconds with the starting time of each interval set to 565 s after the 
starting time of the previous one. 

We htted the light curveof each independent 1130 s time interval in the XMM-Newton 
observations with a constant plus a sine function to map the frequency evolution of the 
mHz QPO accurately. The frequency behaviour of the mHz QPO in the XMM-Newton 
observations is shown in Figure 01 Furthermore, in order to explore the possible mechanism 
responsible for the variation of the QPO frequency, for each XMM-Newton observation we 
divided the intervals with mHz QPO into several segments according to the frequency of 
the QPO (Figured]), and calculated the average frequency and average rms amplitude of 
the QPO in each segment; the details of the segme nts are list e d in T able 121 For the timing 
reduction and analysis of the Obs 4, please refer to iLyn et al.l (l2ni4al ). 


2.2 Spectral data 

For the XMM-Newton observations we applied flag=0 to exclude the events at the edge 
of the CCD and the ones close to a bad pixel to extract the spectra. We produced the re¬ 
sponse matrices using the task rmfgen and the ancillary response hies using the command 
arfgen. For the latter we used the extended point spread function (PSF) model to calculate 
the encircled energy correction. Since the whole CCD was illuminated and contaminated by 
source photons, we extracted the background spectra from another XMM-Newton observa- 


^ We did not exclude the central columns of the CCD to produce these light curves in order to retain enough source photons 
for the detection of the mHz QPO in the power and dynamic power spectra. 

^ We did exclude the central columns of the CCD to produce these light curves in order to measure the frequency and rms 
amplitude of the QPO. 
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Figure 2. Average power spectra of the four XMM-Newton observations of 4U 1636-53 (Obs 1 to 4 from top to bottom) in 
the 0.2 — 5 keV range, calculated from the PN light curve at 1-second resolution. A significant QPO at about 8—10 mHz and 
a weak second harmonic are apparent in each observation. 


tion of the black hole candidate GX 339-4 in timing mode (ObsID 0085680601), based on 


similar sky coordinates and colnmn densit y along the line of sig. 


it. For more detail 


the backgronnd extraction, please refer to 


Hiemstra et ah 


(120 llh and 


Sa.nna et ah 


s abou t 

feoioli . 


We rebinned the source spectra so that we have a minimum of 25 counts per bin, and 
oversampled the energy resolution of the PN detector by a factor of 3. 

We used the heasoft tools version 6.14 for the RXTE data reduction following the 
recipes in the RXTE cook book. We applied the tool saextrct to generate the PGA spec¬ 
tra from Standard-2 data using only events from the best-calibrated PGU2 detector. We 
generated the background using the tool pcabackest and the response files using pcarsp, 
respectively. We extracted the HEXTE spectra for cluster-B events only using the script 
hxtlcurv, and produced the response files using the command hxtrsp. 
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Figure 3. Dynamic power spectra of the XMM-Newton observations of 4U 1636—53. The purple dashed lines indicate the 
time of X-ray bursts in these observations. Each column represents the power spectrum extracted from a time interval of 1130 
seconds with the starting time of each interval set to 565 s after the starting time of the previous one. We oversampled the 
frequency scale by a factor of 100 using the Lomb-Scargle periodogram to display the evolution of the mHz QPOs frequency. 
We fixed the count rate during X-ray bursts and instrument dropouts to the average count rate. The colour bars on the right 
indicate the power at each frequency as defined in the Lomb-Scargle periodogram. 


For Obs 1-3, we fitted the spectrum of each segment from EPIC-PN, PCA and HEXTE 
data simultaneously. For Obs 4, since the RXTE observation covered only segment 1 to 3, we 
combined these three RXTE segments to generate RXTE spectra for segment 4. Besides, we 
noticed that there was a moderately high background flare in the XMM-Newton Obs 3, so 
we applied flare filtering to produce spectra of the total observation. We did not apply flare 
hltering to spectra of the segments in this observation, otherwise we would have excluded 
the segment 3 since most exposure time of this segment was covered by high background 
flare. Alternatively, we only selected PN spectra below 5 keV for spectral analysis to remove 
the possible influence from the flare. 

We first fitted the PN spectra in the 0.8 — 10 keV, and the PCA and HEXTE spectra 
in the range of 10 — 25 and 20 — 100 keV, respectively. We found that there were very 
few source photons in the HEXTE spectra above 30 keV in Obs 3. The parameters in the 
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Table 2. Information for the segments of XMM-Newton observations of 4U 1636—53. Here we give the standard deviation of 
the average frequency since the frequency in segments may change a lot. All errors in the Tables are at the 9 0% confidence 
level unless otherwise indicated. We show the information of the segments in the Obs 4 from iLvu et al.l ll2014ah . 


Observation 

Segment number 

Intervals (s) 

Average frequency (mHz) 

rms amplitude (%) 

Obs 1 

SI 

27120-29380 

10.2±0.3 

1.4±0.2 


S2 

0-2260; 25990-27120; 29380-30510 

9.7±0.1 

1.2±0.1 


S3 

2260-5650 

8.7±0.3 

1.6±0.2 

Obs 2 

SI 

20340-24860 

8.7±0.2 

0.8±0.2 


S2 

16950-20340 

8.2±0.1 

0.7±0.2 


S3 

25990-30510; 31640-32770 

7.6±0.2 

l.liO.l 


S4 

24860-25990; 30510-31640; 32770-35030 

7.1±0.2 

1.2±0.2 

Obs 3 

SI 

565-2825; 8475-9605; 14125-15255 

7.3±0.2 

1.2±0.2 


S2 

2825-6215; 7345-8475; 9605-11865; 12995-14125 

6.8±0.2 

l.liO.l 


S3 

6215-7345; 11865-12995; 15255-16385 

6.1±0.5 

0.9±0.2 

Obs 4 

SI 

0-3954 

8.3±0.6 

0.8±0.1 


S2 

3954-7344 

6.7±0.2 

1.5±0.1 


S3 

7344-11579 

6.9±0.2 

2.3±0.1 


S4 

37310-40800 

9.4±0.3 

1.3±0.1 






Figure 4. Frequency evolution of the mHz QPO in the XMM-Newton observations of 4U 1636—53 (Obs 1 to 2 from left to 
right at the top; Obs 3 to 4 from left to right at the bottom). Each symbol corresponds to the frequency derived from the fit of 
each 1130 s independent light curve using a constant plus a sine function. We divided the QPO cycles in each observation into 
several segments according to the frequency of the QPO. The time intervals within which the QPO showed similar frequency 
were assigned the same segment number as indicated by the rectangles with segment number (Sl-4) in each observation. Since 
the lengths of the segments in Obs 4 are not multiple of 1130 s (see Tablefor details), here we round off the time span of 
each segment in Obs 4 to be a multiple of 1130 s for clarity. 
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fits with HEXTE spectra up to 30 keV were exactly the same as in the hts with HEXTE 
spectra up to 100 keV, and the corresponding reduced chi-square values were much smaller. 
Therefore we decided to use the HEXTE spectra below 30 keV to £t the segments in Obs 3. 

For the reduction of all X-ray bursts detected by the PCA, we produced a spectrum every 
0.25 s during the whole duration of each X-ray burst. We generated the instrument response 
matrix for each spectrum with the standard ftools routine pcarsp, and we corrected each 
spectrum for dead time using the methods supplied by the RXTE team. For each burst we 
extracted the spectrum of the persistent emission just before or after the burst to use as 
background in our fits. 


3 RESULTS 
3.1 Timing results 

In Figure [2] we show the average power spectra of the XMM-Newton observations. We 
detected a strong QPO at about 8—10 mHz and a second harmonic component in each ob¬ 
servation. In Figure |3] it is apparent that the mHz QPO and its second harmonic component 
disappeared right before an X-ray burst in all observations. The QPO in Obs 1 reappeared 
together with its second, and possible third, harmonics ~ 21 ks later after the burst, and 
in Obs 4 the QPO reappeared ~25 ks after the previous burst. In Obs 1 and 3 the QPO 
appeared from the beginning of the observations, while in Obs 2 the QPO appeared 16.4 ks 
after the start of the observation, and then lasted for ~18.6 ks until it became undetectable 
before an X-ray burst. 

In Figure m we show the frequency of the QPO in the intervals of each XMM-Newton 
observation. In Obs 1, the frequency decreased from ~9.9 mHz to ~ 8.5 mHz before the burst, 
it reappeared at a frequency of ~9.6 mHz, then it increased to ~10.4 mHz and returned 
back to 9.6 mHz at the end of the observation. The frequency in Obs 2 started at ~8 mHz, 
it increased to ~9 mHz, and hnally it decreased to below ~7.6 mHz until an X-ray burst 
happened. Interestingly, the frequency in Obs 3 seemed to move randomly around ~6.8 mHz 
before the burst. In Obs 4 the frequency decreased from ~ 8.8 mHz to ~ 6.5 mHz, then 
increased to ~ 7 mHz before the burst, it reappeared at ~9.7 mHz after the burst, and 
then decreased to ~9.1 mHz at the end of the observation. The average frequency and rms 
amplitude of the QPO in each segment are shown in Table |2l we found no clear correlation 
between the frequency and the rms amplitude. 
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Figure 5. Average frequency of the mHz QPO vs. the temperature of the blackbody component in 4U 1636—53. The results 
for Obs 1-4 are marked with black circles, red stars, green squares and blue diamonds, respectively. Here we plot the standard 
deviation of the average frequency in each segment as the error bar of the frequency. The error bar of the temperature is at 
68% confidence level. The black curve in the plot corresponds to the best-fitting power-law model to the data when only the 
error bar of the temperature is used during the fitting process. 


3.2 Spectral results 


We fitted the spectra using XSPEC 12.8.1 (lArnaudlll996h with a model consisting of two ther 


mal components and a Comptonised component (e.g., 


Sanna et ah 


2013 


Lvu et al. 


2ni4bh . 


We used the bbody component to £t the thermal emission from the neutron star surface and 


its boundarv laver. We described tl 

le thermal emission from the accretion disc bv the diskbb 

component ( 

Mitsuda et ah 

1984 

I 

vlakishima et al. 

1986 

. For the Comptonised component 

we selected the model nthcomp 

(Zdziarski et al. 

1996 

Zvcki et al. 

1999 

) to describe the 


non-thermal emission from the inverse Compton scattering process in the c orona, with the 


seed photons coming from the disc (e.g., 


Sauna et ah 


20131: 


Lvu et al. 


2014a 


b|). We used the 


component PHABS to accoun t for the interstellar absorption along the line of si ght, selecting 


the solar abundance table of 


Wilms et al 


fl2000l) and the cross section table of 


Verner et al 


(119961) . We also used the c omponent gauss to account for the iron emission l i ne around 
6—7 ke V in this source (e.g.. 


Pandel et al. 


2008 : 


Ng et al. 


20101 : 


Sauna et al 


20131: 


Lvu et al 


2014bl ). Additionally, we added a multiplicative constant to the model to account for pos¬ 


sible calibration differences between the instruments. Finally, a 0.5% systematic error was 
added to the model. 


For each observation, we assumed the column density, Nj:/, in phabs in each segment to 
be the same and hxed it to be the one in the full observation. We found that the electron 

















































































12 Ming Lyu et al. 

temperature, kTf, in nthcomp, was not well constrained due to the limited exposure time of 
each segment. To overcome this issue we hxed kT^ in the segments to the value derived from 
the hts of the corresponding total observation. We fixed the normalisation of the diskbb 
component in the segments of Obs 2 to the value in the total observation since this parameter 
in the segments could not be well constrained and pegged at its lower boundary. We did 
not use the GAUSS component to fit the segments in Obs 3 since the PN spectra above 5 
keV was removed in order to reduce background flare. The hnal fitting results of the total 
observations and the segments are shown in Table [3] and 01 

We analys ed the spectra of all X-ray bursts in 4U 1636—53 using XSPEC version 12.8.0 
fjArnaudlll996l) . restricting the spectral fits to the energy range 3.0 — 20.0 keV. We fitted the 
time-resolved spectra using a single-temperature blackbody model (bbodyrad in xsPEc) 


times a component that accounts for the interstellar a 
the hydrogen column density, to 0.36x10^^ cm“^ | 


jsorption towards the source and fixing 


Sauna et al. 


2 OI 3 I ). The spectral model 


provides two free parameters: the blackbody colour temperature (Tbb) and the normalization, 
which is proportional to the square of the blackbody radius of the emitting surface, and allows 
us to estimate the bolometric flux during the burst. We then measured the burst fluence 
E^a by integrating the measured fluxes over the burst interval. We used the characteristic 
timescale r = P^bZ/peafe, where fpeak is the peak flux of the burst, to measure the duration 
of the bursts. 

In Figure[5]we show the frequency of the mHz QPO vs. the temperature of the blackbody 
component in the four observations. The frequency covers from ~6 mHz to ~10 mHz, and 
the temperature, kTbb, ranges between ~1.9 keV and ~2.3 keV. We first fitted these data 
with a constant model and got a reduced chi-square of 1.71 {x^/dof=22.23/13). We also 
htted these data with a power-law model, and got a power-law index of 0.16±0.08 (l-a error) 
with a reduced chi-square of 1.50 /dof=18.09/12; for both hts we took into account only 
the errors in the temperature). The F-test probability for these two hts is 0.88, indicating 
that a power-law ht is not signihcant better than a constant. Since Obs 3 was partly ahected 
by solar hares (see §2.2), we also htted the data with a power-law model excluding the points 
coming from this observation. In this case the power-law index is 0.01 ± 0.09. 

In Figure El we show the properties of all X-ray bursts in 4U 1636—53 detected by the 
RXTE satellite as a function of the parameter Sa. Among the ~340 type I X-ray bursts 
detected with RXTE from this sour ce, seven were conhrmed to be associated with mHz 


QPOs: the three bursts reported in 


Altamirano et al. 


(120081) plus the four bursts in the 
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Figure 6. Fluence (upper panel), peak flux (middle panel) and duration (bottom panel) of all X-ray bursts in 4U 1636—53 
detected with the RXTE satellite as a function of Sa (position of the source in the colour-colour diagram; see §2 for more 
details). The black crosses represent all 338 type-I X-ray burs ts, while the blue squares and the red circles represent the seven 
bursts associated with mHz QPOs in lAltamirano et al.l ll200d) and this work. 


XMM-Newton/RXTE observations in this work. The Sa of those 7 bnrsts range from 1.8 to 
2.2, indicating that those bursts appeared around the transitional spectral state. From the 
upper panel of the hgure it is apparent that the fluence of those bursts is relatively high, 
above 40 x 10“® erg cm“^. The middle panel of the figure shows that the peak flux of those 
7 bursts is above 6 x 10“® erg cm“^ s“^, which is much higher than the average peak flux 
of all bursts in this source. In addition, as shown in the bottom panel, those 7 bursts have 
relatively short durations. 


4 DISCUSSION 

We report three new detections of the mHz QPO in the neutron star LMXB 4U 1636-53 
with XMM-Newton. The QPO in all observations in this work disappeared before an type I 
X-ray burst. In the first and the fourth observation, the QPO reappeared at the end of each 
observation after burst, showing a reappearance time of ~ 21 ks and ~25 ks, respectively. 
The QPO in the second observation appeared ~16 ks after the beginning of the observation. 
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Table 3. Fitting results of the observations of 4U 1636—53. Here we give the total unabsorbed flux from 0.5-100 keV in unit 
of erg cm~^ s~^. 


Model Comp 

Parameter 

Obsl 

Obs2 

Obs3 

Obs4 

PHABS 


0.29±0.03 

0.32±0.02 

0.37±0.004 

0.38±0.004 

BBODY 

fcTBB(keV) 

2.25±0.15 

2.03±0.19 

2.01±0.02 

1.96±0.03 


Nor (10-3) 

4.19±0.36 

3.32±0.91 

5.20±0.13 

r 90 + U- tb 
^•^o_0.47 

DISKBB 

kT^bbikeV) 

0.66±0.05 

0.68±0.04 

0.64±0.02 

0.62±0.03 


Nor 

69133 

K9+98 

87l| 

50l|i 

NTHCOMP 

r 

2.52±0.16 

9 qi -(-0.28 

1-1 CO 

o q 
d d 
-h 1 

CO 

CO 

2-49lg:g2 


kTe(keV) 

9-812.9 

q y-(-0.9 
'-0.4 

4.2±0.1 

4 7+2.6 
' -0.4 


Nor 

0.31±0.09 

0.4U_q jg 

0.42±0.005 

0.42±0.03 

GAUSS 

ElineikeV) 

6.4 Ig-^ 

6.55±0.15 

6.75±0.12 

6.80±0.06 


a{keV) 

1 

-*^•^-0.4 

1.5±0.2 

1 o-hO.l 
-^•'^-0.2 

1.21±0.07 


Nor (10-3) 

5-5lkl 

6 

^•^-1.8 

4 q+0-2 

4.1±0.6 


X'^ldof 

202/208 

220/209 

233/208 

242/202 


Total flux (10-^°) 

25.8±0.3 

31.6±0.2 

qq q-(-0.03 

OO. O Q T 

to 

1 + 

o c 

d + 


Table 4. Fitting results of the segments in the observations of 4U 1636—53. Nh and kTe are fixed to the values in the 
corresponding total observations: 0.29x 10^^cm~^ and 9.8 keV for Obs 1, 0.32x 10^^cm~^ and 3.7 keV for Obs 2, 0.37x 10^^cm~^ 
and 4.2 keV for Obs 3 and 0.38xl0^^cm“^ and 4.7 keV for Obs 4, respectively. We also fix the normalisation of the diskbb 
component in the individual segments of Obs 2 to the value in the total observation since this parameter in the segments can 
not be constrained in the fits. For Obs 3, we do not use the GAUSS component at ~6-7 keV since the PN spectra above 5 keV 
was removed to reduce background flare effects. 



BBODY 

DISKBB 


NTHCOMP 


GAUSS 



Obs.Seg 

feT(feeV) 

Nor (10“3) 


Nor 

r 

Nor 

Eli„e(fc<!V) 

crffeeV) 

Nor (10“3) 

/dof 

Ol.Sl 

2.27±0.18 

4.3+0.7 

0.69+0.03 

69+36 

2.52+0.16 

0.30+0.07 

6.65+0.22 

1 0+0-5 

^■^-0.3 

2 9+2-8 
2-0-0.9 

218/209 

01.S2 

2.06+0.16 

4.8 + 1.1 

0.71 + 0.03 

110+^^ 

2.39+0.12 

0.23+0.06 

6 44+0-24 
°-^^-0.04 

1 9+0.6 
1-2-0.3 

2 7+4-3 

2- ^-1.3 

228/209 

01.S3 

2.17+0.14 

4.5+0.6 

0.71 + 0.03 

113+18 

2.39+0.12 

0.23+0.05 

6 43+0-29 
0-4’^_0.03 

1.3+0.4 

3 6+2-0 
cS.O-i 7 

183/208 

02.SI 

2.12 + 0.12 

5.3 + 1.1 

0.74+0.03 

52 

2.42+0.08 

0.38+0.01 

6 57+0-31 

b.57_o 17 

1.2+0.4 

4 5+3-5 
^-^-2.1 

194/209 

02.S2 

2.19+0.25 

2.4+1.0 

0.69+0.03 

52 

2.34+0.06 

0.40+0.01 

54+0.29 

1 r+0-2 

l-®-0.3 

11.2 + 3.1 

222/21Q 

02.S3 

2.04+0.15 

4.7+1.0 

0.74+0.04 

52 

2.39+0.07 

0.39+0.02 

6 7+0.27 
-0.3 

1.4+0.5 

4 9+5-1 
^-^-2.7 

240/210 

02.S4 

2.19+0.12 

5.1 + 1.0 

0.76+0.03 

52 

2.44+0.08 

0.38+0.01 

6 51 +0-29 

o-oi_o.ii 

1.3+0.3 

6.4 + 2.9 

223/210 

03.SI 

2.11 + 0.12 

7.5 + 1.4 

0.65 + 0.03 


2.51+0.10 

0.44+0.07 

- 

- 

- 

140/112 

03.S2 


5.2 + 0.7 

0.65 + 0.03 

134+37 

2.31 + 0.05 

0.36+0.06 

- 

- 

- 

121/111 

03.S3 

1.84 + 0.19 

5.2 + 0.9 

0.66 + 0.04 


2.25 + 0.08 

0.32 + 0.06 

- 

- 

- 

122/111 

04.SI 

'^•'J=_0.07 

4.0+0.7 

0.63+0.03 


2.53+0.09 

0.43+0.07 

6.67+0.12 

1.3+0.2 

6 4+^-0 
°-^-2.0 

234/202 

04.S2 

2.12+0.11 

5.1+0.5 

0.65 + 0.03 

99+24 

yy_4o 

50+0.22 

'^•ou_o,o6 

0.36+0.06 

6.57+0.17 

1 1+0.2 
^'-^-0.1 

4 0+2-3 
^'^ — 1.3 

210/202 

04.S3 

2.03+0.16 

5.2 + 1.0 

0.67+0.04 

118+56 

2.47+0.16 

0.32+0.08 

6 51+0-37 

o-oi_o.ii 

1.5+0.3 

4.7+2.6 

191/208 

04.S4 

^■^^-0.04 

5.4+0.7 

0.65 + 0.02 

= 7+58 

2.50+0.05 

0.40+0.08 

6.74+0.18 

1.1+0.2 

3 7+2-5 

216/202 


and then lasted for ~18.6 ks nntil it disappeared right before a X-ray burst. We found 
no signihcant correlation between the frequency of the mHz QPO and the temperature of 
the blackbody component in the XMM-Newton/RXTE observations. Besides, we also found 
that there was a connection between the mHz QPO and type I X-ray burst properties. 

The mHz QPOs in the four observations in the transitional spectral st ate of 4U 1636-53 


all dis appeared right before X-ray bursts, consistent with the conclusion in 


Altamiranoetah 


(120081 ) that in the transitional state the mHz QPOs disappear only when an X-ray burst 
happens. However, the frequency behaviour before the burst in Obs 3 was different from the 


frequency drift observe d in other observations in the transitional state flAltamirano et al. 


2008 


Lvu et ah! 


2014ah . 
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Altamirano et al. 


(120081) reported a time scale of >15 ks for mHz QPO to reappear after 


the burst when the source was in the transitional state, and a reappearance time of ~ 6 
ks when the source was in the soft state. Considering the two complete QPO reappearance 
times in this work (~ 21 ks in Obs 1 and ~ 25 ks in Obs 4), it appears that the reappearance 
time of mHz QPOs in the transitional state is longer than in the soft state. If that is the 
case, this may indicate that either the X-ray bursts burn more fuel or the fuel accumulation 
process after a burst is slower in the transitional state than in the soft state. 

The frequency of the mHz QPO in Obs 2 hrs t increases and th en decreases for ~15 ks 
until a burst happens. According to the model of iKeek et al.l (120091) . this frequency drift is 
determined by the cooling process of the burning layer, so the frequency drift time of 15 ks 
corresponds to the cooling ti mescale of the bu rning layer. Assuming a burning process with 

({ 2014 ). we found that the thermal cooling timescale 


standard reaction rates from 


Keek et al 


of 15 ks corresponds to a column depth of 2 x 10® g cm“^, which is an order of magnitude 
higher than where the stable H/He burning takes place, therefore, it is likely that the layer is 
heated by the X- ray burst before it cools down to trigger the frequency drift of mHz QPOs. 

According to Heger et al.l (2007) (see their eq. 11), and assuming that the thickness of 
the fuel layer remains more or less constant during the mHz QPOs cycle, the neutron-star 
temperature (that we identify with the blackbody temperature in our hts) should be propor¬ 
tional to (z//rh)“^, where v is the frequency of the mHz QPO, and m is the mass accretion 
rate onto the neutron star. Assuming that the total unabsorbed flux (F) is proportional to 
mass accretion rate, we htted kTss vs. v/F with a power law relation and got a power law 
index of 0.11±0.05 (l-cr error) with a reduced chi-square of 1.4 ( y® /dof=17.33/12). This 
number is signihcantl y different from th e index of —2 predicted by 

The simulations of! 


Keek et al. 


Heger et ahl (120071) 


(I 2 OO 9 I ) suggest that there should be a correlation between 


the frequency of the mHz QPO and the temperature of the blackbody component: As the 
burning layer gradually cooles down, the frequency of the oscillations decreases by tens of 
percents until an X-ray burst happens. To investigate how our observational constraints 
on the power-law index comp are to mult i -zone simulations, we considered the hydrogen- 


accreting models presented by 


Keek et ahl 


(120141) . We determine the power-law index [3 for 


the different series of models where marg i nally stable burning was present by using the 


oscillation period reported by 


Keek et al. 


(120141 ) and the time-averaged effective surface 


temperature of the models. For most models /3 ~ 0.13, whereas the simulations with a 
decreased rate of the (a, 7) ^®Ne reaction yield (3 = 0.22. A study of pure-helium burning 
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models yields {3 ~ 0.04 flKeek et al. 


2009 


for M = O.SMEdd, without rotationally induced 


mixing). Our observational constraint on (3 is consistent with these values. Interestingly, the 
predicted change in (3 for the reduced (a, 7 ) ^®Ne rate is similar to the la uncertainty in 
our measurement. Future improved observations of mHz QPOs hold a promise to constrain 
the uncertain nuclear reaction rates of CNO breakout reactions (e.g., iDavids et al.ll 201 ll) 


processes ( 

Fisker et al. 

2006; 

CooDer & Naravan 

2006; 

Fisker et al. 

2007; 

Parikh et al. 

2008; 

Davids et al. 

2011 

). Besides nuclear reaction rates and accretion composition, the stability 


been exhaustive ly studied, such 


transition depends on a series of parameters that have no 
as the base luminosity a nd rotationally in duced mixing flKeek et al.ll2009l) . as well as the 


eff ective surface gr a vity 


Hegeretaf 


2003). 


Galloway et al.l (120081 ) showed that the long bursts with slow rise and decay were likely 


mixed H /He rich bursts, while f ast b ursts with shorter rise 


Helium. 


Sugimoto et al. 


f l984ll 


and 


Galloway et al. 


Ime were likely to burn primarily 


(120061 1 investigated the distribution of 


peak fluxes of bursts in 4U 1636—53; these results show that some faint photospheric radius 
expansion (PRE) bursts have relatively low peak fluxes. According to their interpretation, 
during the expansion of the photosphere, a H-rich outer layer is ejected, exposing the un¬ 
derlying helium layer below. For the normal PRE bursts, the bursts reach the FEdd,He after 
the ejection of the outer H-rich layer, while the faint PRE bursts have insufficient energy 
to drive off the H-rich layer, reaching only pEdd.H instead of pEdd.He- In this work, all seven 
PRE bursts associated with mHz QPOs are short, bright and energetic, suggesting that they 
reach pEdd.He and hence belong to He-rich bursts. Therefore, it is very likely that there is 
potential connection between the mHz QPOs and the He-rich PRE bursts when the source 
is in the transitional state. 
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